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method, 简称 为 DG 方法 ) 是 由 Reed 和 Hill 在 解决 中 子 输 运 问 。 度 的 DG 方法 。 近 年 来 ,CFD 工作 者 的 不 懈 努 力 丰 富 和 发 
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摘 要 : 为 了 求解 三 维 欧 拉 方程 ， 对 隐 式 时 间 离 散 格 式 间 断 有 限 元 方法 进行 了 研究 。 根 据 间 断 Galerkin 有 限 元 方法 思 
想 ， 构 造 内 迭代 SOR-LU-SGS 隐 式 时 间 离 散 格 式 ， 结 合 当地 时 间 步 长 技术 、 多 重 网 格 方法 ， 实 现 了 三 维 流 场 的 计算 。 
数值 计算 了 ONERAM6 机 可 、 大 攻 角 类 前 缘 三 角 翼 以 及 DLR-F4 翼 身 组 合体 的 亚 声 速 绕 流 问 题 。 结 果 表明 ,加 入 SOR 
内 和 迭代 步 的 LU-SGS 隐 式 算法 具有 较 大 的 优势 ， 相 较 于 GMRES 算法 所 占用 的 内 存 少 且 收敛 速度 相当 ， 是 LU-SGS 算 
法 的 3 倍 以 上 。 针 对 三 维 算 例 , 具有 较 好 的 稳定 性 和 较 高 的 收 竹 速度 ,能够 给 出 准确 的 流 场 信息 。 与 原 方 法 相 比 ，SOR- 
LU-SGS 方法 无 论 是 在 迭代 步 数 上 还 是 在 CPU 时 间 上 ， 效 率 均 有 明显 提高 ， 适 合 于 三 维 复 杂 流 场 计算 。 
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Implicit discontinuous finite element algorithm for solving three-dimensional Euler equations 


Duan Zhijian, Zhang Tong, Qin Mengmei, Ma Xinrong 
(School of Mathematics & Information Science, Xianyang Normal University, Xianyang Shaanxi 712000, China) 


Abstract: In order to solve three-dimensional Euler equations, this work investigated the implicit time discretization 
discontinuous finite element algorithm. According to discontinuous Galerkin finite element method ideology, this paper 
established SOR-LU-SGS time discretization discontinuous Galerkin finite element method for three-dimensional inviscid flow 
fluid coupled with implicit time discretization scheme, local time step technique and multi-grid method. Then, it calculated the 
classical three-dimensional examples, including subsonic flow around the ONERA M6 wing, delta wing with tip, and DLR-F4 
wing-body geometry. The results show that the SOR-LU-SGS implicit algorithm has more advantages than GMRES with less 
memory, and the convergence rate is more than triple times of the traditional LU-SGS algorithm. For the three-dimensional 
examples, SOR-LU-SGS scheme not only has good stability and high convergence rate, but also can give accurate flow field 
information. Compared with the original method, the present scheme improves the efficiency both in iteration steps and in the 
CPU time, which is very suitable for solving three-dimensional complex flow field. 
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DG 方法 和 FVM 混合 算法 求解 了 三 维 N-S 方程 , 得 到 了 DG 方 


引言 
只 法 在 处 理 高 超声 速 流动 时 具有 优势 。 奉 凡 等 人 [9 针对 三 维 
间断 有 限 元 方法 (discontinuous Galerkin finite element RANS 方程 研究 了 叶轮 机 机 械 内 流 问 题 ， 发 展 了 一 种 高 效 高 精 


题 时 而 提出 的 。 由 于 DG 方法 具有 在 网 格 单元 内 部 可 以 构造 任 ”DG 方法 ， 并 且 广 泛 应 用 于 各 类 工程 问题 0 。 


意 高 阶 的 精度 ， 易 于 处 理 边界 和 实施 并 行 计算 等 诸多 优点 ， 所 由 于 工程 问题 的 复杂 性 , 提高 计算 效率 是 CFD 工作 
以 被 认为 是 CFD 领域 最 具 前 景 的 算法 之 一 P31。 近 年 来 在 航空 、 的 目标 。 提 高 计算 效率 行 之 有 效 的 方式 主要 包括 : a. 采 用 隐 


展 了 


生 创 


向 为 


格式 
加 速 


海洋 、 气 象 、 石 油 勘探 等 诸多 领域 得 到 了 广泛 应 用 。1989 年 以 ” 间 推 进 格式 或 结合 加 速 技术 的 显 式 时 间 推 进 格式 ; b. 并 行 计算 。 

来 ，Cockburn 等 人 [7 发 表 了 一 系列 研究 成 果 ， 构 建 了 求解 双 并行 计 算是 发 展 的 必然 趋势 ， 本 文 主要 针对 第 一 点 进行 
守恒 律 方程 的 框架 。Bassi 等 人 8 利用 混合 有 限 元 方法 ， 推 常用 的 显 式 时 间 推 进 格式 为 显 式 Runge-Kutta 多 步 

六 到 了 层 流 N-S 方程 的 求解 问题 。 痪 立新 等 人 0% 成 功 地 利用 “0949， 辅 之 以 多 重 网 格 方法 、 当 地 时 间 步 长 、 残 值 光 顺 等 
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录用 稿 


技术 模拟 无 香 流 动 取得 了 巨大 的 成 功 。 但 是 当 处 理 大 雷诺 数 等 
刚性 问题 时 ， 计 算 效 率 迅 速 下 降 。 因 此 ， 采 用 仅 依 赖 于 物理 流 
动机 理 的 隐 式 时 间 离 散 格 式 成 为 更 好 的 选择 ， 其 中 以 LU-SGS 
格式 0 和 GMRES 格式 09 最 具 代 表 性 。 
本 文 改进 了 LU-SGS 隐 式 时 间 离 散 格式 ， 构 造 了 适合 于 
维 双 曲 守恒 律 方程 的 隐 式 DG 方法 。 通 过 大 攻 角 尖 前 缘 三 角 恤 
和 DLR-F4 器 身 组 合体 的 亚 声速 绕 流 问 题 对 算法 进行 了 验证 ， 
计算 结果 与 实验 结果 吻合 良好 。 通 过 ONERA M6 机 翼 跨 声速 
绕 流 问题 对 GMRES、LU-SGS 格式 和 Runge-Kutta 多 步 格式 及 
改进 算法 进行 了 分 析 比 较 ， 改 进 算法 的 效率 得 到 了 明显 提升 ， 
且 节 约 了 计算 成 本 。 


二 


ll 


1 ”三 维 Euler 方程 
不 考虑 体积 力 和 外 部 热源 ， 可 压缩 Euler 方程 为 


OQ(x,?) 
gr VOR) (1) 


: Q 、 分别 为 守恒 变量 和 通 量 函 数 ， 分 别 定义 如 下 : 


这 里 : P,e,P 分 别 为 气体 密度 、 单 位 体积 总 内 能 及 压强 ; 


Ui 为 x 方向 的 速度 0; 为 Kronecker 符号 。 加 入 理想 气体 状态 


方程 : 
p=(7— pe 1g 
(r=D( puuj),h=et~, 


形成 封闭 方程 组 ， 其 中 y=1.4 为 空气 比 热 比 。 


间断 有 限 元 方法 的 构造 


基于 DG 方法 定义 函数 的 有 限 元 空间 如 下 : 


vr={v el :vaelVl oeQ), © 


其 中 : VY 为 局 部 函数 空间 , 是 p(p =b2,…) 次 多 项 式 集合 ， 信 


为 子 域 空间 
假设 求解 域 中 的 近似 解 为 2, ， 在 方程 式 (2) 两 侧 乘 以 试探 
ee 


5], GOdQo- 上 人 F.(0) .2 ‘dO+ 


。 G3) 
上 F,(Q,)-n®,dT =0, 0 e a 


其 中 ，@@, 为 试探 函数 中 的 近似 解 ，=(n., 太 ,7.) 为 单元 边界 
,的 单位 外 法 向 量 ， 并 且 Q, 、 中 , 满足 
Q; =2,00070%), D, = ,8,070%), (4) 


其 中 ， 多 (Go 是 己 阶 基 函 数 ， 且 同时 满足 式 (3)。 
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pa (9) 
| FQ) -ng Var =0, 1< jsN 
] 黎 曼 通 量 函 数 玉 (QO7,QX,n) (Qi 、QX 为 单元 界面 两 侧 音 
元 解 ) 代 替 五 (Q@,)-n; 项 ， 联 立 式 (3)~(5) 可 得 如 下 方程 组 : 
MQ = Res(O). (0 
dt 
为 了 提高 计算 效率 和 节约 存储 空间 ， 本 文采 用 标准 正 交 晶 
函数 简化 矩阵 结构 。 式 (6) 为 关于 时 间 ; 的 常 微分 半 离散 方程 组 ， 
相对 于 显 式 推进 格式 ， 隐 式 推进 格式 往往 是 无 条 件 稳定 的 ， 医 
此 本 文采 用 隐 式 推进 格式 。 


3 ” 隐 式 时 间 推 进 格 式 


相对 于 显 式 推进 格式 而 言 ， 隐 式 推进 格式 计算 量 小 、 稳 定 
性 好 ,是 目前 CFD 常用 的 计算 方法 。 相对 于 多 次 迭代 ,传统 的 
LU-SGS 格式 仅 实施 了 一 次 对 称 Gauss-Seidel 和 迭代， 在 一 个 时 
间 步 中 计算 效率 明显 下 降 ， 整 体 计算 效率 较 低 。 为 此 ， 本 文 类 
合 内 迭代 SOR 方法 ,联合 多 次 对 称 技术 , 得 到 了 SOR-LU-SGS 
隐 式 时 间 离 散 格式 。 

将 式 (6) 运 用 差分 离散 及 线 化 处 理 后 ， 得 到 NN 维 线性 方程 
组 MAO" = Res”" ，NN 为 总 网 格 数 ， 此 时 M 严格 对 角 占 优 。 分 
解 方程 组 MAO"” = Res" 可 得 

(L+D+U)AQ”" = Res” (7) 
其 中 : 工 、D 、D 满足 式 (8)。 
4 1 = es, j<i 


N 
p=- T+ > 和 “Hi 十 


i 
At /a 


1 


y 


rs, 5 (8) 


i 


7 -| 加 


y 


00; 
略 去 高 阶 小 量 ， 化 简 (7) 式 可 得 
(D+DD (D+U)AQO" = Res”. 


1| oF (|,)) Pk 
3 rs, j>i 


进一步 采用 两 步 推进 : 
(D+ DAQ' = Res”, 
(D+U)AQ”" = DAQO”™. 
加 入 内 和 迭代 SOR 格式 后 ， 具 体 实现 过 程 如 下 : 
AQ =(1-0) A + oD 四 -2 (ON , "Pls AQ, 
(9) 


zi rl 1 网 -中 oo 
1 


AC0 =(1-o)AC + oD" 四 -> 站 


7 2 
oF (©,) 


1 (9) 
-到 全 | 


这 里 ， 内 和 迭代 次 数 由 式 (11) 决 定 


af(2) 
60， 


n, -|% | |AS, AQ; 


(10) 
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录用 入 
|soe -ace of <s, al) e=6.5x105， 总 网 格 量 为 290 万 。 
取 &=0.1， 松 弛 因子 6。 满足 0 < w < 2。 
最 终 和 欠 代 收敛 后 ， 可 得 
2” =2 +AQ" 
此 外 ， 文 中 的 通 量 函 数 采 用 RoeP9 格 式 ， 并 且 配 合 当 地 时 


间 步 长 技术 、 多 重 网 格 技术 进行 加 速 处 理 。 


4 数值 实验 
4.1 算法 验证 

算 例 1: 尖 前 缘 三 角 细 的 大 攻 角 亚 声速 流动 

计算 状态 : 马赫 数 Ma = 0.4 , 攻 角 =18 , 总 网 格 量 为 500 
万 。 


mn Experiment results 
Sotfware(CFX) results 
SOR-LU-SGS results 


pe 


| 2 
9 


8 
RE 


ee CK lj 
图 2 沿 展 向 剖面 压力 系数 曲线 


图 1 给 出 了 三 角 翼 表面 网 格 。 图 
结果 以 及 实验 结果 在 不 同位 置 处 的 比较 。 可 
计算 的 结果 与 本 文 算法 
维 


采用 

的 流动 为 非 定常 流动 ， 因 此 在 0.8 之 后 两 者 具有 差异 。 
算 例 2: DLR-F4 辟 身 组 合体 的 亚 声速 流动 
计算 状态 : 马赫 数 Ma = 


2 是 本 文 算法 、CFX 计算 
以 看 出 ， 利 ) 
的 结果 保持 一 致 ， 基 本 与 实验 值 
双 曲 守恒 律 方程 求解 的 涡 破裂 位 置 在 0.75 附近 , 此 后 


0.6 ， 攻 角 w =0' ， 实 验 雷 
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DLR-F4 器 身 组 合体 表面 网 格 


(a) 18.5% 


position=0.844 


Experimen 


(b) 84.4% 


1 results 


机 辟 不 同 


] CFX 
切合 。 


虹 向 剖面 的 压力 系数 出 


5 DLR-F4 辟 身 组 合体 表面 


j 压 力 系数 云图 


图 3 


给 出 了 DLR-F4 机 避 表 面 网 格 ， 


区 


4 给 出 了 机 械 


展 向 


201804.01446v1 


chinaXiv 


录用 稿 


不 同 站 位 (18.5% 和 84.4%) 压 力 系数 与 实验 值 对 比 结果 。 可 以 看 
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SGS 方法 相 比 ， 提 高 了 3 倍 以 上 ， 与 Runge-Kutta 格式 相 比 提 


出 计算 结果 与 实验 值 基本 吻合 。 然 而 计算 值 稍 高 于 实验 值 ， 这 
主要 是 因为 机 辟 安 装 角 较 大 以 及 采用 了 超 临 界 翼 型 ， 从 而 导致 
压力 分 布 对 经 型 厚度 很 敏感 。 图 5 给 出 了 全 机 表面 压力 系数 分 
布 。 
4.2 ” 算 例 求解 与 分 析 

算 例 1: ONERA M6 机 辟 跨 声速 流动 

计算 状态 : 马赫 数 Ma =0.84， 攻 角 =3.06 ,实验 雷诺 数 
Re =11.7x105 。 


图 6 ONERA M6 表面 网 格 


Computational results 


(a) 20% 


position=0.95 


Experiment results 


Computational results 


(b) 95% 
图 7 ONERA M6 机 辟 不 同 展 向 位 置 的 压力 系数 曲线 


图 6 给 出 了 ONERA M6 机 辟 表 面 网 格 ， 图 7 给 出 了 不 同 
展 向 位 置 处 的 压力 系数 对 比 曲线 。 可 以 看 出 与 实验 值 基本 吻合 。 
图 8 中 可 以 看 出 本 文 的 计算 结果 很 好 的 捕捉 到 Y 型 激 波 结构 。 

图 9 和 10 中 分 别 给 出 了 四 种 时 间 离 散 格式 收敛 历史 。 显 然 本 
文 构造 的 SOR-LU-SGS 方法 具有 较 好 的 收敛 速度 , 与 传统 LU- 


高 了 近 6 倍 ， 基 本 接近 于 GMERS 格式 ， 但 本 文 算法 所 需要 的 
存储 量 却 远 低 于 GMRES 格式 。 


0.463 


0.033 


-0.396 


-0.826 


图 8 ONERA M6 机 辟 表 面 压力 系数 


A CT TT 


ee Runge-Kutia method 
------- SOR-LU-SGS method 
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图 9 不 同 算法 时 间 收 敛 历史 


eee Runge-Kutta method 
= SOR-LU-SGS method 
号 MHRES method 
LU-SGS method 


Res(p) 
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图 10 ”不同 算法 和 狗 代 步 数 收敛 历史 


隐 式 格式 一 般 需 要 存储 系数 矩阵 ， 导 致 存储 量 较 大 。 本 文 
加 入 内 迭代 LU-SGS 算法 相对 于 显 式 Runge-Kutta 格式 只 需要 
增加 一 个 Yx N 的 数组 来 存储 Jacobi 矩阵 和 对 角 和 矩阵 ， 即 计算 
量 和 存储 量 分 别 为 COV212)+O(CV212) 和 O(M2) 。GMRES 隐 
式 时 间 推 进 格 式 的 计算 量 和 存储 量 分 别 为 
O(N)+O(N”)+Olm) 和 OOzV)+OCV2)+OCOzD) ， 随 着 m 和 
N 的 增 大 ，GMRES 格式 的 存储 量 远 高 于 本 文 算法 的 存储 量 。 


5 结束语 
本 文采 用 SOR-LU-SGS 隐 式 时 间 离 散 格 式 , 结合 当地 时 间 
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步 长 技术 和 多 重 网 格 方法 ， 构 造 了 求解 三 维 双 曲 守恒 律 方程 的 
隐 式 间断 有 限 元 方法 。 由 计算 结果 可 以 得 出 ， 本 文 算法 具有 较 
高 的 收敛 速度 和 较 好 的 稳定 性 ， 针 对 于 ONERA M6 机 辟 跨 声 
速 绕 流 流 场 计算 ， 收 敛 速度 相对 于 传统 的 LU-SGS 格式 提高 了 
3 倍 以 上 , 相对 于 显 式 Runge-Kutta We 6 倍 , 效率 基 
本 接近 于 GMRES 格式 ， 且 所 需 存储 量 较 低 。 总 之 ， 本 文 方法 
适合 于 三 维 复杂 流 场 计算 ， 也 为 进一步 求解 可 压缩 Navier- 
Stokes 方程 打下 了 良好 的 基础 。 
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